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Csl ; Abstract 

jrt ■ We study the dynamics of a tree-level A-type atoms driven by a coherent train of short, non- 

overlapping laser pulses. We derive analytical non-perturbative expressions for density matrix 

O ' 

p^ \ by approximating pulses by delta-function. We demonstrate that depending on train parameters 






several scenarios of system dynamics are realized. We show the possibility of driving Raman 
transitions between the two ground states of A-system avoiding populating excited state by using 
the pulses with effective area equal to 27r. The number of 27r-pulses needed to transfer the entire 



d ' population from one ground state to another depends on the ratio between the Rabi frequencies of 

O ' two allowed transitions. In the case of equal Rabi frequencies, the system can be transferred from 

^~>, one ground state to another with a single 27r pulse. When the total pulse area differs from 2tt and 

Qh| the two-photon resonance condition is fulfilled, the system evolves into a "dark" state and becomes 

ly-. ■ transparent to subsequent pulses. The third possible scenario is the quasi-steady-state regime when 

tH- . neither the total single pulse area is equal to 27r nor the two-photon resonance condition is fullfilled. 

m , 

^^ ' In this regime the radiative-decay-induced drop in the population following a given pulse is fully 

p' 

-Y-\ ' restored by the subsequent pulse. We derive analytical expression for the density matrix in the 

O 

^ ■ quasi-steady-state regime. We analyze the dependence of the post-pulse excited state population 

in the quasi-steady-state regime on the train parameters. We find the optimal values for train 

^ ' parameters corresponding to the maximimum of the excited state population. The maximum of 

H ■ 

C^ ' the excited state population in the steady state regime is reached at the effective single pulse area 

equal to vr and is equal to 2/3 in the limiting case when its radiative lifetime is much shorter then 

the pulse repetition period. 

PACS numbers: 32.80.Qk, 42.50.Hz 



I. INTRODUCTION 



The frequency combs (FC) generated by the trains of ultrashort laser pulses [l| have been 
actively developed over the past 10 years. Recerrtly a flber-laser-based FC with 10 W average 

power was demonstrated [2] with the prospects of further scaling of technology up to 10 kW 
average power. The spectral coverage of combs has been extended from optical to the ultra- 
violet and mid-IR spectral range |3|-l5| . High resolution quantum control via the combination 
of pulse shaping and frequency comb was shown in 6|-[9|. The experiments on line- by- 



line addressing have been done 



7|. These rapid technological developments enable novel 



applications 



computing [12 



in precise metrology [l|, llO| , atomic and molecular spectroscopy [ll| , quantum 



13| 



manipulating external and internal degrees of freedom of atomic and 



molecular systems 14|4l6|. In this paper we explore the dynamics of three-level A- type atoms 
Fi g. P a)) interactin g w ith coherent train of ultrashort laser pulses. The dynamics of two- 



17H20|, three 



-|l6. 



21 



Ci 



22| and multi-level systems J14l] driven by such trains has been actively 



investigated over the past decade. In particular there were proposals for Doppler cooling of 
atoms based on two-photon transitions driven by ultrafast pulse trains |23|, optical pumping 
and vibrational cooling of molecules by femtosecond-shaped pulses [l5| and rotational cooling 
of molecules by chirped laser pulses |16 |. 
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FIG. 1: (color online ) Energy levels of A-system and positions of frequency comb teeth. The 
comb is Doppler shifted in the atomic frame moving with velocity v. 



The analytical expression for the density matrix of a two-level system interacting with 



the pulse train was obtained in 20|, l22| However, in many cases the atom can not be ap- 
proximated as a two-level system because the excited state decay to intermediate sublevels. 
As a practicle example of A system, we consider the ground states of group III atoms. Their 
ground states are composed of two fine-structure sublevels nPi/2 and nP^i2 and the decay 
to the intermediate level can not be neglected |24| . 

Recently there was a series of works studying the dynamics of three-level atoms interacting 



with a train of ultrashort pulses |2ll. l25l |26 | 



In particular, accumulative effects in the coherence of three-level atoms excited by 
femtosecond-laser frequency combs was studied in 2l|. There authors obtained pertur- 



bative iterative solution 
trapping was studied in 



or the density matrix of three-level system. Coherent population 



25 



26| . Perturbative analytical iterative solution for the density 



matrix in a weak field limit was presented. 

Here we derive the analytical non-perturbative expressions for the density matrix of a 
A-system driven by coherent trains of ultrashort laser pulses. Our work can be considered as 
extension of earlier works for a two-level system driven by the pulse train 20|, |22| . As in our 
previous work on two-level system [20|, we use the model of delta- function shaped pulses. 
Using the derived equations we study dependence of system dynamics on the parameters of 
the pulse train. 

For the pulse-train-driven A-system there are two major qualitative effects: "memory" 
and "pathway-interference" effects. Both effects play important role in understanding of 
multilevel-system dynamics driven by the pulse train. 

The system retains the memory of the preceding pulse as long as the population of the 
excited state does not decay between subsequent pulses. Then the quantum-mechanical 
amplitudes driven by successive pulses interfere and the spectral response of the system 
reflects the underlying frequency- comb structure of the pulse train. If we fix the atomic 
lifetime and increase the period between the pulses, the interference pattern is expected to 
"wash out" , with a complete loss of memory in the limit of large decay rates. This memory 
effect is qualitatively identical to the case of the two-level system, explored in Ref. 20|. 

The "pathway-interference" effect is unique for multilevel systems. The excited-state 
amplitude arises from simultaneous excitations of the two ground states. The two excitation 
pathways interfere. The "pathway-interference" effect is perhaps most dramatic in the CPT 
where the "dark" superposition of the ground states conspires to interfere 



regime 



25 



destructively, so that there is no population transfer to the excited state at all. 

We show that in a particular case when the integer number of FC teeth fits into the 
energy gap between the two ground states and the total single pulse area differs from 27r 
the system evolves into a "dark" superposition of ground states and becomes transparent 
to the following pulses. The ratio of the populations of two ground states in this regime is 
determined by the ratio of corresponding pulse areas. This effect is commonly referred to 



as a coherent population trapping (CPT) |25l428| 



We also show that when the total single pulse area (defined as the geometric sum of 
individual pulse areas 9i, 62, corresponding to two different transitions, = \/6j + 6^ ) is 
equal to 27r then Raman transitions can be driven between the two ground states avoiding 
the excited state. The number of pulses needed for complete population transfer from one 
ground state to another depends on the ratio between two pulse areas 61, 62- In a particular 
case of equal pulse areas the entire population can be transferred from one ground state to 
another by a single pulse. 

Finally, we derive analytical expression for the density matrix of a system in a steady 
state regime realized for finite decay rate of the excited state and the total single pulse 
area 7^ 27r. In this regime the radiative-decay-induced drop in the population following 
a given pulse is fully restored by the subsequent pulse. We analyze the dependence of the 
quasi-steady-state post-pulse excited state population on the FC parameters. 

This paper is organized as follows. In section I we derive general non-perturbatuve recur- 
rent equation for density matrix of A-system interacting with a coherent train of ultrashort 
laser pulses. In section II we enumerate main parameters characterising interaction of A- 
system with a pulse train. In section III we study different scenarios of the system dynamics, 
each realized for certain combination of parameters. Finally, conclusions are drawn in Sec. 

m 

II. ANALYTICAL SOLUTION OF THE OPTICAL BLOCH EQUATIONS FOR A 
DELTA-FUNCTION PULSE TRAIN 

In a typical setup, a train of phase-coherent pulses is generated by multiple reflections of 
a single pulse injected into an optical cavity. A short pulse is outcoupled every roundtrip 
of the wavepacket inside the cavity, determining a repetition time T between subsequent 



pulses. At a fixed spatial coordinate, the electric field of the train may be parameterized as 

E(t)=£Ep^cos(a;,t + $^)(?(t-mT), (1) 

where e is the polarization vector, E^ is the field amplitude, and $m is the phase shift. The 
frequency lOc is the carrier frequency of the laser field and g(t) is the shape of the pulses. 
We normalize g{t) so that max \g{i)\ = 1, then Ep has the meaning of the peak amplitude. 
While typically pulses have identical shapes and $,m = ^{rriT), one may want to install an 
active optical element at the output of the cavity that could vary the phase and the shape 
of the pulses. 

We are interested in a dynamics of three-level A-system, interacting with the train ([1]), 
see FigHJ A-system is composed of the excited state |e) and the ground states \gi) \g2) 
separated by A12; the transition frequencies between the excited and each of the ground 
states are cjggj, Ueg2 correspondingly. The optical Bloch equations (OBE) for the relevant 
density matrix elements (populations Pee, Pgigi, Pg2g2 ^^^ coherences peg and Pge, j = 1,2) 
read 



-IP: 



ee 



-^/m[fieg,pegj , (2) 



i=i 



2 Q* 

Pegj = -l^Peg, + i Y^ -^ iPeeSjp - pg,g^ ) , (3) 

99,9,1 = ^jj'ljPee + T^i^lg Pg,e - ^eg.Peg^,)- (4) 



yjyj' JJ iji c^ ' r) 

The time- and space-dependent Rabi frequency is 



N-l 

C 



^eg, {z, t) = fif "' Y,g{t+-- ,riT)e-*('^^(*)-^^*-*™) , (5) 

rra=0 

where 5j = Uc — ^egp ^c = ^c/c and z is the atomic coordinate). The peak Rabi- frequency 
VlF^"" = -^(e|D ■ ^l^f) is expressed in terms of the dipole matrix element. Eqs. ([21 |3l H]) were 
derived using the rotating wave approximation. Notice that the energy gap between the two 
ground states can be expressed in terms of individual detunings: A12 = 62 — 61. 

Notice that as long as the duration of the pulse is much shorter than the excited state 
lifetime and the repetition time, the atomic system behaves as if it were a subject to a per- 
turbation by a series of delta-function-like pulses. In this limit, the only relevant parameter 



affecting tlie quantum-meclianical time evolution is tlie effective area of tlie pulse 

oo 

e, = nf''' I g{t)dt, (6) 

— CO 

and i7^^" g(t) — )■ 9j5{t) in all the previous expressions. As an illustration, we may consider 
a Gaussian-shaped pulse, g{t) = e^* /^^p. In the limit r^ <^ T, this pulse is equivalent to a 
delta-function pulse lim e~* '^"^p — > y/2nTp6{t), as both pulses have the very same effective 



area 9j. 



Now we turn to finding the solution of the OBEs for a coherent train of delta-function 
pulses, 

Af-l 2 

Pee = -IPee " ^ ^(^ " ^^) E^^^'^^ [e-*(^=^W-^^*"*-)pe,,] , (7) 

m=0 jf=l 

. N-l 2 

Pe., = -|pe,, + \Y.^it- tj Y. ^pe^^^-W-^^*-*-)(p,e<5,p - P,,,J, (8) 

m=0 p=l 

. AT-l 

%.,, = '^.y7,Pee + I E ^(^ - ^T)(^,,e^(^-W-^.'*-*-)p,^,e " ^,e-'(^-W-^^*-*™)pe,^, ). (9) 

m=0 



We will distinguish between pre-pulse (left) and post-pulse (right) elements of the density 
matrix, e.g., (p™); and (p™)^ are the values of the excited state population just before and 
just after the {m + ly^ pulse. Below we relate these values at each pulse and between the 
pulses. Starting from given initial values of p and applying a recurrent procedure we may 
find p at later times. 

Delta-function pulses cause abrupt changes in density matrix elements at points tm, = mT. 
Between the pulses, however, the dynamics is simple as it is determined by the spontaneous 
decay. 

This leads to the following time evolution between the pulses {mT < t < (m + 1)T)) 

Peeit) = {p:Xe-'\ (10) 

Pe,,(t) = (pS,)/-"*, (11) 

P,^,At) = (p-^,)^ + 6,,,^ (pTX (1 - e-^*). (12) 

Further, we may neglect the spontaneous decay during the pulse, since for a typical 
femtosecond pulse 'jTp ^ 1. Then the OBEs in time interval (mT — e < t < mT + e), 

6 



e — )■ O"*", may be recast in the form p = —iS{t — mT) [am, p], where [a^, p] is a commutator 
and the matrix a^ reads: 



/ 



ar 



^le^'^iW ^2e^''2(*) 
0ie-'''iW 

02e-^''2(*) 



(13) 



Here 



r]i{t) = kcZ-Sit-<^{t). 



(14) 



The matrix notation corresponds to the following enumeration scheme for matrix elements 
of p 



Pee Pegi Peg2 



P 



Pgie Pgigi Pgig2 
\Pg2e Pg2gi Pgig-iJ 
The exact analytical solution of this equation is p (t) 



(15) 



U\t)p{mT - e)U{t) = 
, with T being the time- 



W{t){p"')^U{t), where U{t) = Texp ia^ y^j,_J{t' - mT)dt' 

ordering operator. Thus the pre- and post-pulse elements of the density matrix are related 

by 



(P" 



Am \P )i A^ 



(16) 



where 



A. 



' cos ® 

-i sin ® sin ;i^e"*''^*^*™^ 



— zsin ^ sin xe*''^ *•*"'■* 
cos^ X + cos y sin^ x 



— z cos X sin y e*^2 ^*™ ^ * 
sin^ f sin(2x)e"*'^i2im 



\ — zcosxsm ^e 



e^-ir,2(t„) _ g|j^2 e sin(2x)e*^i2t™ 



sm 



e 



(17) 



X + cosfcos^X J 



Here 6 = a/^i + 6*2 is the total single pulse area defined as the geometric sum of individual 
single pulse areas of the two transitions and x = arctan ( |^ j determines their ratio. At this 
point, by combining Eq. ( jTOj) and Eq. (TT6|) one may find time evolution of the density matrix 
over a single repetition period; apparently, by stacking these single-pulse and free-evolution 
propagators, one may evolve a given initial p over duration of the entire train. In Fig. [2] 
we show results of such calculation for the excited state population of a heavy atom (atom 
remains at rest). 

In Fig. |2]the atom is initially in the lowest ground state \gi). The lifetime of the excited 
atomic state is 15 ns, and the decay rates are equal: 71 = 72 = 7/2 . If the frequency gap 



AnT=1200.l9-^2K 




FIG. 2: Time-evolution of the excited state population in A-system interacting with a coherent train 
of laser pulses. The carrier frequency is resonant with the transition between the lowest ground 
and the excited state. The energy gap between the two ground states is A12 = 27r x 300 GHz and 
the decay rates are 71 = 72 = 27r x lOMHz and B\ = Bi- (a) Comparison of two-photon-resonant 
(dark purple line - A12T = 1200 x 27r) and off-resonant (pink line - A12T = 1200.19 x 27r) regimes. 
The effective single pulse area is = ^, ^i = 02- (b) The effect of the pulse area for a fixed value 
of A12T = 1200.19 X 27r: dark purple curve corresponds to G = vr and the solution for = ^ is 
shown in pink. 

between the two ground states A12 is commensurate with the pulse repetition rate (27r/T) 
(see Fig. |2] (a), dark purple line), then the system evolves into a "dark" superposition of 
ground states and becomes transparent to the pulses. 

III. CHARACTERISTIC DIMENSIONLESS PARAMETERS 

To streamline the analysis we introduce dimensionless parameters, characterizing pulse- 
train cooling of the A-system. 

(i) The ratio of the pulse repetition period and the lifetime of the excited state 



/i = 7T, 



(18) 



This parameter will in particular characterize the spectral profile of the post-pulse 
excited state population. 

(ii) Single-pulse areas Oj for the two transitions \gj) — )> |e), j = 1, 2. We will also employ 
two related auxiliary parameters: the angle determining the ratio between the single 
pulse areas 9i, 62 

X = arctan(0i/02) (19) 

and the effective single-pulse area 



Q = ^el + el (20) 

(iii) Branching ratios, based on the decay rates of the excited state to the two ground states 

&i = 7i/7, fc2=72/7- (21) 

Certainly 61 + 62 = 1- 
(iv) Number of teeth fitting in the energy gap h/\i2 between the two ground states 

K = Aia/wrep • (22) 

Notice that k generally is not an integer number. When it is integer, the two-photon 
resonance conditions are satisfied and the system evolves into the dark state. 

(v) Doppler shifted phase ^\M offsets between subsequent pulses defined as 

r)^=7]i{t)-7]i{t + T) = {k,V + 5i)T + (j), Ti2 = m + 2TTK. (23) 

Here v is the atomic velocity and is the carrier-envelope phase offset between subse- 
quent pulses, i.e., = $m+i — ^m in Eq- dl])- These phase parameters will be used to 
characterize the spectral profile of the excited state population. As shown below the 
density matrix of a system is a periodic function oifii,fi2. The two phases are always 
related as 

% —fji = 27r/t = 2nAi2/u!rep ■ 



(vi) Residual detunings Sj, j = 1,2, between \gj) levels and the nearest FC modes in 
the reference frame moving with the atom. In general, 6i = rj^/T + l-nnxjT and 
82 = fi2/T + 27Tn2/T, where integers rij = 0, ±1.. are chosen to renormalize the residual 
detunings to the interval —Urepf^ < 5j < Wrep/2. 

IV. SYSTEM DYNAMICS 

Below we show that the system dynamics is mostly determined by four parameters n, x, 
G, 7. Depending on these parameters the following four scenarios may be realized. These 
different regimes are covered in individual subsections of this section. 

(a) Dark state (CPT) regime is realized for finite decay rate 7, when the integer number 
of FC teeth fits into the energy gap between the two ground states (/t = 0, 1..). Here 
the system evolves into a stationary superposition of two ground states ("dark" state), 
which is transparent to the pulse train. 

(b) Stimulated Raman transitions between the two ground states (avoiding populating 
the excited state) are observed in the A-system when the effective single pulse area 
is O = 27rn, n = 0, 1.. and the decay of the excited state within the pulse can be 
neglected {'-/Tp ^1). If initially the system is in one of the ground states, then the 
excited state remains unpopulated after each new pulse and the system evolves as 
a time-dependent superposition of two ground states \gj). Pulses lead to, discussed 
below, abrupt change of coefficients in this superposition. As shown below, at some 
special choice of x? the entire population can be transferred from one ground state to 
another by a single = 27r pulse. The decay of the excited state can be neglcted for 
the number of pulses estimated as A^ ~ l/7Tp- 



(c) If the lifetime of the excited state is much longer than the pulse repetition period T 
then for a number of pulses, N <ti l/{'~fT), the dissipation can be neglected. In this 
case, if the effective single pulse area is not a multiple of 27r, (0 7^ 27rn), the population 
in A-system oscillates between all three states. This is the transient regime preceding 
the quasi-steady-state regime. 

(d) The quasi-steady-state regime (QSS). After N ^ l/(7rp) ^ 1/(7T) pulses the system 
evolves into a saturated regime. In this regime, the same fraction of population is 

10 



driven to the excited state by each pulse, so the maximum value of (Pee)^ is reached 
at the moment of time just after each pulse. Between the pulses the excited state 
population exponentially decays to the ground states and reaches its minium value 
just before the next pulse. These minimum and maximum values of the excited state 
population do not depend on the sequential number of the pulse. 

A. "Dark" state (CPT) 

When the energy gap between the two ground states is commensurate with the distance 
between modes in a FC (k = 0, 1, .), the two-photon resonance condition is fulfilled 26|, |27|, 
and (similar to the case of two CW sources) the Hamiltonian posesses stationary "dark" 
state. Here the atom is in a superposition of two ground states, described in the interaction 
picture by the wave function 



/dark 



cos(x)l^i) -sin(x)|^2). (24) 



Once in the stationary state, the system dwells in it unless perturbed (e.g., pulse train 
parameters change). As a result, the system becomes transparent to the pulse train. This can 
be also explained by the distructive interference between quantum probability amplitudes 
of the transitions \gj) ^ |e) at k = Au/uJrep for the system in a "dark" state. 

The fact that the superposition flMl) is an eigenstate can be observed from the fact that 
the density matrix, corresponding to ([2i 



dark 



^0 ^ 

sin^x -'-^ 

sin(2x) 2 



(25) 



^ T^ cos X ) 

commutes with the time-evolution operator Am flT7|) . Notice that the dark state f l24|) does 
not depend on the branching ratios &i,&2- The two-photon resonance (k = 0,1..) is a 
prerequisite for the existence of a stationary state in A-system. Below we show that the 
"dark" state can be avoided for a large number of pulses (A^ ~ l/{pfTp)) if the effective pulse 
area is = 2?™, n = 1, 2. 
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B. Stimulated Raman transitions between the two ground states 

When the effective single pulse area is = 27rn, n = 1.., and the decay of the excited 
state during the pulse can be neglected {-yTp <^ 1), the system oscillates between the two 
ground states, avoiding populating the excited state altogether. If the radiative decay during 
each pulse can be neglected and the excited state population {pl^)^. is zero, then analytical 
expression for the time-evolution operator after the A^-th pulse can be obtained as the 
product AAr..A3A2Ai, where the operator Am is defined by Eq. (fTTj) . Knowing the time- 
evolution operator, one can express the wave function (which initially was in the lowest 
ground state) after the A^-th pulse as 



,N 



where 



ndsp - Cg^^N\gi) + Cg2,N\92), (26) 



C,uN = -^YTTT-, {{-ifFMv) + FM-V)) , (27) 

sm (2x) 

-i{N-l)TTK 

Cg,,N = -—^^ ii-l^F.^i^) - n^{-^)) , (28) 



FaNi^) = e^^- ,^u.^-.u.,nr.j....^j _^ ^^Q) 



1 + csc2(2x) (cos(y9 — cot(7rfi;) sin;/)) 
j^^ (cos v? — cot(7rK) sin Lp) 



F,^{^) = e^^^ ^ ^ ^ ' ^^ 2, (30) 

1 + csc^(2x) (cos(y9 — cot(7r/T;) simp) 

sinyj = sin(7r/T;) cos(2x) cosv9>0. (31) 

In particular, when n = 1/2 + 2n {n = 0, 1..), one has: 
C.,. = e- (i±<^ co= {n (I - 2,)) - ,i^ =in (iV (f - 2, 



(32) 



The system which is initially in one ground state can be transfered to another ground 
state hj N = 2k pulses when x = ii;(2(^ ^ ^ 1) ^^^ hj N = 2k — 1 pulses when 
X = 4(2fc-i) (2(^ ^ ~ 1)5 ^ = 1,2..., / = 0,..k — 1. In a special case of equal pulse 
areas (for example, 9i = 62 = V2-K, {x = ^/4)) the entire population can be transferred 
from one ground state to another by a single O = 27r pulse. If initially the excited state 

12 



was populated, then pee either remains constant if there is no decay to the lower states or 
becomes distributed between the oscillating populations of the two ground states if there is 
a decay of excited state to any of the ground states. 

It is worth highlighting the difference in meaning of the 27r-pulse in two- and three-level 
systems. In a two-level system the 2tt pulse would drive the entire population to the excited 
state and then return to the ground state by the same pulse simultaneously. In the case 
of three-level system one could explain vanishing excited state population at the end of the 
27r pulse (if it was zero before the pulse) in a similar fashion the same pulse drives the 
population to the upper state and then back to the superposition of the two ground states. 
The nature of this pjpcess is different from the well-known Stimulated Raman Adiabatic 
Passage (STIRAP) |1J], involving two CW sources with slow- varying amplitudes and equal 
detunings between carrier frequencies and transition frequencies. In our pulsed laser case 
driving the population between the two ground states avoiding the excited state is not 
affected by the difference in detunings 6i, 62- In the limiting case when the excited state is 
metastable 7 — > 00, the conclusions made here can be generalized for slow varying-envelope 
pulses as long as the conditions for 6 and x remain fulfilled. 

C. Transient regime 

During initial sequence of A^ ^ 1/(7T) pulses the decay of the excited state can be 
neglected and the density matrix evolves as 

Ptransient = AjV-i.. Ai AqPoAqA^.. A^, (33) 

where po is the initial density matrix. In this regime, the A-system oscillates between all 
three states. At k = 0, 1.. the wave function describing the system after the A^-th pulse (if 
initially all the population is in the ground state \gi)) can be expressed as 

\'^)Zansient = CgM + Cg^\g2) + Ce|e), (34) 

Ce = -«smxsm-— , (35j 

Cg^ = cos^(x) + sin^(x) cos -^, (36) 

Cg, = -sin(2x)sin^^. (37) 



(38) 
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During the transient regime the "dark" state is not reached yet even if the two-photon 
resonance condition is fulfilled |25| . 

D. Quasi-steady-state regime 



Similar to the case of two kicked coupled damped pendula [29|], the system eventually 
reaches saturated regime, where the radiative-decay-induced drop in the population following 
a given pulse is fully restored by the subsequent pulse. We will refer to this behavior as the 
quasi-steady-state (QSS) regime. As shown in the Appendix, the QSS regime allows for a 
fully analytical solution. Since pee (t) = pee it + riT), in the QSS regime, pre- and post-pulse 
values (p^); ^ do not depend on the pulse number m and we denote these values as (Pee); r- 
Furthermore, because of the radiative decay, {ple)i = ^""'"^ {ple)r- "^^^ general solution 
for the density matrix in the QSS regime can be obtained from a system of linear algebraic 
equations derived from condition: (p^) , = (p^^^) ,, where (p^) is the re-normalized density 
matrix (see the Appendix for details). The solution is fully analytical, however it is unwieldy 
and here we present its simplified form obtained for equal pulse areas 6i = ^2 (x = 7r/4, = 

Solution for the arbitrary pulse area is given in Appendix. 
The post-pulse value is 

■yT 

(pa = ^sin2Msm^|, (39) 



where 



D ^ ' 2 



D = [hi cos rj^ + 62 cos r^g) I 4 cos — — sin^ 77 ~ ^ cos(27rfi;) cos'^ — 
—2 (62 cosfj^ + hi cosr/2) I sin"^ — + cos^ — 
-|-2sin(27r/t)(62 sinr/2 + hi sinr/^) cos^ - — 



-2 cosh — I sin^ — + cos — sin^ (ttk) ) . (40) 



The equation flA9p is symmetric with respect to the swap of the ground state labels. The 
dependence on the phase offset r/^ (after substituting % = ^1 + ^t^k) is the result of inter- 
ference between the elementary responses of a system to subsequent pulses (the persistent 
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"memory" of the system). Particularly, when 7T — )■ 00, the excited state completely decays 
between the pulses and the interference factor vanishes (the "memory" is erased), 



4sin^ (tik) 



~^ tan^ ^ + '^^"'('^7) ' 
At equal branching ratios bi = b2 = 1/2 the equation (1A8I) can be simplified further 

e 2 sm (ttk) sm ^ 



(p: 



4L)' 
D' = ( cos (ttk) cos {f]i + ttk) ( cos^ (ttk) cos"^ ( T ) ~ '^'^^ ~o ) ~'~ 

7T^ ^_4 /0^ , 2 /0\ . 2 



cosh (^^ J (^sin^ (^- J + cos^ \^-j sin^ (vr/.) J J . (42) 

We start analyzing the expression (H2]) for the QSS value of the excited state population 
by studying its spectral profile. In Fig. |3]^a) we plot the dependence of the excited state 
population (pl^,)^ on the phase offset r/^^. As an illustration we choose the following set of 
parameters: k = 0.12 + 2n {n = 0, 1..), 61 = 62 = 1/2, 7T = 1/4, O = 7r/3 (for example, 
7 = 27r X 10 MHz, T = 3.9804ns, A12 = 27r x SOOGif^). 

The periodic structure mimics the FC spectrum. The maxima of {ple)r ^^^ reached at 
r/]^ = —mod^TTK, 27r) + 2nn, n = 0, ±1.., that is at symmetrical values of detunings (5i) = 
- p2)^pj = -mod{7TK,27T)/T. 

The optimal ratio of the residual detunings (61/62)^ ^ can be defined as the value corre- 
sponding to the maximum post-pulse excited state population. At equal pulse areas 61 = 62 
and equal branching ratios 61 = 62 = 1/2 this optimal value is equal to —1. Generally, it 
depends on the ratio between branching ratios and the ratio between individual pulse areas 
61/62. In case when one of the frequencies Uegj is nearly resonant with the nearby FC tooth 
and another is not resonant with any of the FC modes, one would observe accumulation of 
the population in the state which is "less coupled" (not resonant). This is the reason why 
at equal pulse areas (^1 = 62) the accumulation of population in one of the ground states is 
reduced by setting the detunings 5i = —82 = —mod{TCK,2Tc)/T. The detunings have to be 
of opposite signs to avoid the two-photon resonance which would drive the system into the 
"dark" state. 

At different pulse areas and equal branching ratios 61 = 62 = 1/2 the decay rates of the 
excited state to both ground states are equal but the rate of the repumping of population 
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FIG. 3: (a) The dependence of the QSS excited state population on the phase offset parameter 
rji at fixed values of parameters 7T = 1/4, B = vr/S, 61 = 1/2, 62 = 1/2 and k = 0.12. (b) 
The dependence of the optimal ratio of the residual detunings \{6i/62)opt\ on the parameter X) 
corresponding to the maximum value of the post-pulse excited state population. The solid purple 
curve is obtained for equal branching ratios (61 = 62)- The optimal ratio |(^i/^2)opt| for the case 
when the branching ratios are varied with the parameter x as 61 = sin^ x = ^i/©^! ^2 = cos^ x = 
^2/0^ is shown with the dashed pink curve. 

from a certain ground state \gj) to the excited state depends on the pulse area 9j and the 
residual detuning 6j. When 61 = —62, one would expect that the state corresponding to 
smaller pulse area 6j accumulates the population and the excited-state population vanishes. 
To mitigate this effect, the state of smaller effective pulse area 6j has to be closer to the 
resonance with FC tooth (smaller detuning 6j) than another one, corresponding to the 
larger effective pulse area. Therefore, one would expect that at different pulse areas and 
equal branching ratios, 61 = 62 = 1/2, the optimal ratio between the detunings {61/62)^ ^ at 
which the post-pulse excited state population has its maximum value, grows with increase 
of the parameter tanx = 61/62- It is worth noticing that at A12 ^ Ueg-^ equal branching 
ratios mean equal dipole matrix elements, entering the definition of Rabi frequency ([5]). In 
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order to obtain different pulse areas 61^2 at equal branching ratios 61 =62, one would need 
additional pulse shaping, to modulate intensities of different FC teeth. 

When the pulse areas scale proportionally to the square roots of branching ratios (all 
the teeth have the same intensity) and 61 7^ 62? the situation is different. In this case, 
the ground state of smaller pulse area is "less coupled" to the excited state, but the decay 
rate of the excited state to this ground state is slower. We found that in this case the 
optimal ratio {81/52) decreases when increasing the ratio tanx = |^. In Fig. [3]^b) we plot 
the dependence of the ratio {81/52) ^ ^ as a function of the parameter x at different values 
of branching ratios &i,&2- The solid purple curve is obtained for equal branching ratios, 
61 = 62 = 1/2. The dashed pink curve was drawn assuming that the branching ratios vary 
with the parameter x as 61 = sin^ x, 62 = cos^ x- 




FIG. 4: The dependence of the quasi-steady-state value of post-pulse excited state population 
on parameter k = ^yi/iOrep at h\ = 62, 1^ = 1/4, Q = vr/3. Panel (a) rji = 2Trn. Panel (b) 



Vi 



-7r/10±27rn,n = 0,1... 



Next we study the dependence of the post-pulse excited state population value on the 
parameter k, i.e., the ratio between the ground states energy gap and the pulse repetition 
frequency Urep- In Fig. H] (a,b) we plot the dependence of the excited state population (Pee)r 
( 142|) on the parameter k = Ai2/ujrep at different values of r/j^. Both curves exhibit periodic 
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pattern which mimics the periodic spectrum of the pulse train. The dips at integer values 
of K correspond to the CPT regime with zero excited state population. Manipulating the 
pulse repetition rate Urep stretches the positions of FC modes in the frequency domain and 
consequently the residual detunings 5j, between the frequencies Ueg^ and nearest teeth. 

The maxima of population in Fig. H] (a,b) remain the same when increasing the value of 
^12/ ^rep (that is increasing the value of T). This can be explained by the fact that at fixed 
value of the parameter 7T, the value of the excited state population (pgg)^ depends on 5j 
and ujrep only through the ratios 5j/urep and Au/cOrep, see Eq.( H2|) . 

Notice that the profile in FigJD^b) is asymmetric, while the one on FigJH (a) is symmetric 
with respect to the integer values of k. To explain this assymetry we parameterize f/i and 
K as f/i = 27ma + St], k, ^ nt, + 6k,, where free parameters t] and k, are constrained as 
^ 6ri < 2tt and ^ 6k < 1. When the ria-th harmonic is resonant with the frequency u^g-^ , 
61] = 0. Then different values of k = nh±6K correspond to the frequency Ueg2 being red(blue) 
detuned with respect to the {ua — nb)-th mode by 6KUrep- Corresponding values of residual 
detunning ^2 is 62 = ±6KijJrep if 6k < 1/2 and =f(1 — 6K)ijJrep if 6k > 1/2. Flipping the sign 
of 62 (at 61 = 0) does not affect time-evolution of the system, causing the dependence of 
iPee)r '^^ ^ (Fig- 13(a)) at f^i = to be symmetrical with respect to /t = ra, n = 0, 1... 

If rj^ differs from the integer multiple of 27r, fj^ = 27ma + 6t], (where 6t] < 2-7?), then 
the Ua-th FC harmonic is detuned from the frequency Wegi by -^cOrep- At k = nh± 6k both 
frequencies Ueg^ generally do not match any of the FC modes. Different values oi k = na±6K 
correspond to different detunings ^2 at fixed value of 6i, causing the dependence (Fig. Hl^b)) 
of {plf,)r o^ K at f]i = — 7r/10 to be asymmetric. However the "translational" symmetry with 
respect to the shift k = K±n, n = 0,1.. still remains. 

As we showed for fixed k the optimal value of residual detuning, is (5i) = 
—mod{TTK,2TT)/T. Now we would like to vary k in order to optimize Pee further. One 
can find that this optimal value oi k = k°^^ can be expressed as 

K°P* = larccos(x), (43) 

where a; is a root of the following algebraic equation: 



„4„_40 oo^„„.T.7T.:^4 , i«_ o^2/.___0 



16x^ COS* - - 32a; cosh -- sin* - + 16 cos 2a;^ 4 cos - + 3 cos 9 + 9 =0. (44) 

4 2 4 2 V 2 \ \ I 

18 



At fixed values of the decay rate 7, pulse area 9, and the frequency gap between the two 
ground states A 12, the equation (jUj) is a self-consistent equation for T. 



E. Maximum post-pulse excited state population in the quasi-steady-state regime 

In previous subsection we found that the maximum of the post-pulse excited state pop- 
ulation (Pee)r ^^ reached at optimal residual detunings 61 = —62 = —mod{n°^^ /2, 1)/T and 
optimal parameter k = k°p* determined by Eq. (1431) . 

Now we would like to vary the pulse area O to optimize this maximum. In Fig. \5\ we plot 
the dependence of {ple)r '^^ ^^^ effective single pulse area 9. Different curves correspond to 
different values of parameter /i = 7T. The values of (Peg)^ were calculated at the optimal 
value of K, determined by Eq. (jUj) for each and p = 7T. 




FIG. 5: The dependence of the quasi-steady-state values of the post-pulse excited state population 
(Pee)r ^^ effective single pulse area at different values of /i = 7T: /i = 10 (dashed pink line), 
/i = 1/2 (dashd blue line), fi = 1/100 (solid purple line) and optimal parameters 771 = —k°^^/2, 
where k°p* is obtained from Eq. (j43p . 

From Fig. |5]we see that the maximum values of the excited state population are attained 
at 6 = TT -|- 2?™. Substituting 6 = vr in the equation (Hij) . one finds: 

1 



^opt 



+ n, 



n 



0,1... 



(45) 



For these values of k and G the excited state population and the fractional momentum kick 
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are: 



{plMQ = vr, r/i = -TTK, K=^) = ^e^^/V cosh(7r/2) (46) 

(47) 

The spectral resolution of the excited state population vanishes as O — )■ tt and k = Kopt- 
The maximum of {ple)r i^ three-level A-system (with 61 = 62 = 1/2, 61 = 62 = \2,t{) 

is reached at 7T ^1. Its value is 2/3 that is different from the case of two- level system, 

where the maximum excited state population in the quasi-steady state regime is 1, as it was 

shown in our previous work 201]. 

In case of unequal pulse areas 9i 7^ 62 and branching ratios 61 7^ 62, at O = vr and 7T ^ 1, 

the three-level A-system, which is initially in the ground state \gi)-, eventually reaches the 

QSS with post-pulse excited state population expressed as 

2sin^ 7rfi:sin^(2x) 
(61 - 62) cos(2x) + 1 + 2 sin^ ttk sin^(2x) ' 



{Ple)r (0 = ^, iT » 1) = 77 ..„_,..,. , . , .:. 2 r^T^TJ- (48) 



.2^, _ el , _ _2 



If the branching ratios vary as 61 = sin X = '^1^2 = cos x = '^1 the Eq. ( I48l) does not 
depend on the value of x, 

ipieX (0 = ^> IT » 1, h = sin^ x) = ,^"fZ^ - (49) 

2 sm TTK + 1 

The maximum value of {ple)r i^ ^^^^ ^^^^ i^ reached at k = 1/2 + n, n = 0, 1.. and is equal 
to 2/3. In case if 61 = 0, the system which starts in the ground state \gi) obviously stays 
unperturbed. In case ii 61 = n, x = ^/2, 61 = 1 the system is reduced to a pair of coupled 
levels \gi) and |e). Then, the maximum population inversion and fractional momentum kick 
are equal to 1. 

V. CONCLUSION 

In this paper we studied the dynamics of a three-level A-type system driven by a train 
of ultra-short laser pulses. General analytic expressions for time-evolution of the density 
matrix were obtained. Several regimes of system dynamics can be realized depending on the 
train parameters. 

In particular, when the two-photon resonance condition l^yil'^re-p = 0, 1.. is fulfilled, the 
system evolves into a stationary "dark" state where it becomes transparent to the pulses. 
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In the limiting case when the total pulse area is a multiple of 27r, the "dark" state 
is avoided. In this case, regardless of the pulse repetition rate and the decay rate of the 
excited state the post-pulse excited state population vanishes. The system oscillates between 
the two ground states avoiding populating the excited state alltogether. In a special case of 
equal pulse areas the entire population can be transfered from one ground state to another 
by a single 6 = 27r pulse. 

At finite excited state decay rates, the system eventually reaches the quasi-steady-state 
regime which is similar to the saturated regime in a system of two kicked coupled damped 
pendula. In the QSS regime the radiative-decay-induced drop in the population following a 
given pulse is fully restored by the subsequent pulse. 

We derived analytical expression for the density matrix in the QSS regime, neglecting the 
decay during the pulse. The post-pulse excited state population has a periodic dependence 
on the Doppler shifted phase offset between the subsequent pulses. This periodic pattern 
reflects the frequency comb spectrum and strongly depends on the ratio between the pulse 
repetition period and the excited state lifetime. 

In a particular case when the pulse repetition period is much longer then the excited state 
lifetime, the interference between subsequent pulses vanishes and the spectral dependence 
of the excited state population mimics the spectral profile of an individual pulse. 

In the opposite case when the excited state lifetime is much longer then the pulse repeti- 
tion period and the single pulse area is small O ~ 7T the pulse train acts on a system as a 
collection of narrow-band CW lasers with individual frequencies corresponding to different 
FC modes. 

At a given pulse area the maximum of excited state population is reached at some optimal 
ratio of residual detunings between the frequencies of the two allowed transitions and the 
nearest FC teeth. This optimal value depends on the effective single pulse area, branching 
ratios and the ratio of individual pulse areas. At equal branching ratios and equal pulse 
areas the optimal residual detunings have the same absolute value and opposite sign. The 
single pulse area corresponding to the maximum population inversion is equal to vr. In 
case when the ratio of individual pulse areas are determined by the ratio of corresponding 
dipole matrix elements only, the absolute maximum of the QSS population inversion in the 
saturation regime (reached at = tt and 7T ^ 1) does not depend on the ratio of these 
dipole matrix elements and is equal to 2/3. In this case the optimal residual detunings are 
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(5i 2 = ±'^^, This result is different from the case of two-level system, where the maximum 
population inversion in the saturation regime was equal to 1. 

Appendix A: Density matrix in the saturation regime 

Here we derive the value of the density-matrixreached in the saturation (quasi-steady- 
state) regime. The pre- and post-pulse elements of the density matrix at the A^*"^ pulse are 
related by Eq. lTTB]) 



Introducing the unitary transformation 



(Al) 



ip)r = U^NiP)Lr^N, 



where 



UN 



I .■ '7l(tjv)+'?2(*iv) 

' e 2 








„- '?2(tiv)-')l(tjv) 

e 2 














' 
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tjv)-12 


tjs, 




e* 




2 




/ 



one can rewrite (lAlll as: 



(A2) 



(A3) 



(p^)^ = Ao (p^), At. (A4) 

The quasi-steady-state density matrix (p^) = (p*)^ can be obtained from the system of 
linear equations 

(p").= (p"-^).- (A5) 

The general equation for the post-pulse excited state population can be expressed then as 



(p: 



4e^^ sin^ (I) sin2 (vr^:) sinh^ (2^ ) sin2(2x) 



D 



(A6) 



where 



D = Det ({Di, D2, Ds, D4, D,, Dq, Dj, D^}) 



(A7) 
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Di = |1— e' — (1 + cos x) sin — -, — sin — - cos2x, 0, 0, sin —sin 2%, 
sin 9 sin x, sin 6 cos x, 0} 

rri /~\ /~\ 

D2 = {26ie^~ sinh h 4sin^ xsin^ — (1 — sin^ — (1 + cos^x)), 

2<S> 2 f 2^ \ 2© f 202 

—4 sin -rsin x I sin — cos2x + 1 J , 0, 0, — 2sin — sin2x I 1 — 2sin -r ^^^ ^ 

e / 2© 2 \ 2© 

—2 sin — sin X I 1 — 2 sin -r sin x J, 2 sin -7 sin — sinxsin2x, 0|, 
D3 = {0, 0,1 — e^" cos f]i — 2 sin^ — cos^ x? sin^ — sin 2x, 0, e^" sin f]i, 0, sin — cos x}, 

2© 221/_ N 2©2 ^ / - 

D4 = {0,0, sin — sin 2x, 1 — 6 2 cos(?7i + 27r/t) — 2 sin -r sin X)0,0,e2 sin(r7i + 27r/T;), 

-sm — smxl, 

D^ = {sin^ — sin2x(3cos^ sin^ — cos2x), — sin"^ — sin4x, 0, 0, 

4© . 4© © /^ . 2 2©~ 

— sm — cos 4x + cos - — cos ztck, — cos x sm — - I 1 — 4 sm x sm -r ; , 

© f. 2 2© \ /NT 

sm— -smxl4cos xsm - — 1 I , sm(27rK)|, 

©/ 2 / 2 X ©N © /^ 2© \ iT _ 

Dq = {sin X sin —(cos x " (1 + cos x) cos — ), sin — sinx I 2sin — cos2x + 1 1 , — e 2 sin?7i,0, 

•© /^1.•2©•2^ 2T 62 ^-2 

sm — cos X I 1 — 4 sm —sm X I 5 ^^ ^ cos rji + cos — cos x + cos k3 sm X; 

e\ sin2x „. 
cos H — cos — I , 0|, 



6,2©/ 2 X N © /^ 2© A 

Dj = {2 cos X sin —(sin — (1 + cos x) ~ l); ~ cosxsin — I 1 — 2sin — cos2x),0, 

— e^~ sin(f/i + 27rK), sin — sin x I 1 — 4 sin^ — cos^ X ) ? ~ sin^ — sin 2x I 2 cos hi 

.jr ,_....( © ^\ . 2 



— e 2 cos(r7i + 27rK) + I cos — — cos 9 j sin x + cos 6, 0}, 
Dg = {0, 0, — cos X sin —, sin — sinx, — sin(27rfi;), 0, 0, cos cos(27r/i;)} 



In some limiting cases the general equation (1A6P can be simplified further. 

(a) At ex = e2- 
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where 



^Pee)r = "^ sill^ {txk) sill^ ^ , (A8) 



D'^ = (61 cos r]i + 62 cos 772) I 4 cos — — sin^ 77 "~ 2 cos(27r/i:) cos^ — 
—2 (62 cosr];^ + 61 cos 772) I sin'' — + cos^ — - 
+2sin(27r/T;)(62 sinfjg + &i sinr/j^) cos^ 

7T r 4 , 2 0.2 



—2 cosh — I sin — + cos — sin (ttk) J . (A9) 



(b) At 61 = sin^ X, 62 = cos^ x- 



(Pee)r = 8e~ sin^ — sin^ TT/t/D* 

D = 8 cos 2x j 4 sin"^ — + sin^ — cos 271 k J sin ttk sin (f]i + ttk) + 

+ cosTT/tcos (r^i + TT/t) I 4 COS — (cos27r/t — 5) + (cos© + 3)(3cos27r/t + l) — 
— 16sin^ — sin^ 7r/i:cos(4x) J — 
-4 cosh ( — ) ( 4 cos^ — cos 27rfi; + 2 cos — -cos© -5 j . (AlO) 



(c) At © = tt: 



rri 

iPee)r ~ 16e^~ sin^ vr/t sinh(-— ) sin^ 2x/ D^ 

Dc = sinh ( — j (cosr^i (2(2(6i - 62 + 2)cos27r/t + 61 + 562Cos47r/€+ 1) + 

+ cos2x(2(86i - l)cos27r/t + 6i - 1562Cos47r/t + 1) + 8sin^ (ttk) cos4x(fei - 2 - 
—862 cos27r/T;) + 4sin^ 7r/i:cos6x(&2 cos27r/i: — 61)) — 8sinf]i (4^2 sin 47rfi; sin^ x+ 
+ sin27r/T:cos^ x(^2(cos4x + 3) — 4cos2x))) + 
+4 sinh(7r) (2(6i - 62) cos 2^-2 cos 27r/t sin^ 2x - cos 4% + 3) (All) 
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(d) At 9 = TT, and bi = sin^ X; ^2 = cos^ x- 



D'^ = cosf/i(8cos27rfi; + 4cos2x — cos(4x) + 5) — 16(cos27r/t — 2) cosh I -— I + 
2 sinful sin2'7rK(4cos2x + cos(4x) — 1) + 8 sin^ x cos(f/i + Atik). (A12) 
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